LOCAL HOLDER CONTINUITY FOR SOME DOUBLY 
NONLINEAR PARABOLIC EQUATIONS IN MEASURE 

SPACES 



EURICA HENRIQUES AND ROJBIN LALEOGLU 

Abstract. We establish the local Holder continuity for the nonneg- 
ative weak solutions of certain doubly nonlinear parabolic equations 
possessing a singularity in the time derivative part and a degeneracy in 
the principal part by using the method of intrinsic scaling. We work 
bX) in measure spaces equipped with a doubling non-trivial Borel measure 

^ ' supporting a Poincare inequality. 



1. Introduction 

We consider the regularity question for the nonnegative weak solutions of 
the doubly nonlinear equation 

d(u q ) 

-4r^ - V • (|Vd p - 2 Vu) = 0, p>2 and < q < 1. (1.1) 
at 

This particular equation and similar ones appear in the modeling of turbu- 
lent filtration of a gas or a liquid through porous media and non-Newtonian 
fluids (see [9]). In the case q = 1, the equation becomes the heat equation 
for p = 2 and the degenerate p-Laplace equation for p ^ 2. 

In addition to the double nonlinearity equation (jl.ip possesses, it has 
a degeneracy coming from the fact that its modulus of ellipticity |Vit| p-2 
vanishes at points where |Vu| = 0, and a singularity due to the term u 9-1 , 
enclosed in the time derivative part, which blows up at points where u = 0. 
By using the method of intrinsic scaling, (see [8j [26] ) , we show that weak 
solutions are locally Holder continuous in the setting of a measure space 
with a doubling non-trivial Borel measure supporting a Poincare inequality. 
For other successful uses of intrinsic geometry, see also [U El El Q31 [20] . 

There exists a dichotomy in the behaviour of weak solutions of (jl.ip . 
For large scales, when the infimum is considerably smaller than the oscil- 
lation, the singularity in the time derivative part dominates. However, for 
small scales the equation becomes essentially the singular p-Laplace equa- 
tion. Both cases are treated via the method of intrinsic scaling but each 
requires the construction of a suitable intrinsic geometry on its own and the 
different choice of certain test functions. 

This work can be regarded as a natural continuation of [21] (see also 
[19] ) which considers the case q = p — 1, corresponding the so-called the 
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degenerate Trudinger's equation. Unlike Trudinger's equation (see [25| [T7]). 
in our case the Harnack's inequality is not available. Therefore, despite 
the similarities with the Holder continuity analysis for this equation the 
techniques employed are quite different. 

The same regularity issues were considered by Porzio and Vespri, Vespri 
and Ivanov, see [22], [271 DSL f° r equations that are equivalent to (jl.ip pro- 
vided the solution is strictly away from zero (see also |10t [Tl] ) . 

The rest of the paper is organized as follows. In section 2 we provide the 
basic definitions and facts on metric spaces that are needed for our goals 
and present also the main result. Section 3 contains the construction of 
the iteration argument and fundamental estimates, namely the energy and 
logarithmic estimates; sections 4 includes the continuity analysis. The case 
corresponding to the parabolic behaviour is easier and will be treated at the 
end of the 4th section. 



2. Preliminary material and main result 

Let fi be a Borel measure and U be an open set in M. d . The Sobolev space 
H 1,p (U\n) is defined to be the completion of C°°{U) with respect to the 
Sobolev norm 

l/p 



\ u \\l,p,U 



j (\u\ p + \Vu\ p )dfj\ 



A function u belongs to the local Sobolev space H^ P (U;fi) if it belongs 
to H 1,p (^l;n) for every compactly contained subset Q, of U. Moreover, the 
Sobolev space with zero boundary values H Q ' P (U ; (j,) is defined as the com- 
pletion of C^°(U) with respect to the Sobolev norm. For more properties of 
Sobolev spaces, see e.g. [13] . 

The parabolic Sobolev space LP(t 1 ,t 2 ;H 1 ' p (U;fi)), with t x < t 2 , is the 
space of functions u(x, t) such that, for almost every t € (t\, t 2 ) the function 
u(-,t) belongs to H 1,p (U;/jl) and 



I 2 [ (\u\ p + \Vu\ p )dv < oo, 
Jti Ju 



where we denote du = dfidt. 

The space L p oc (t\,t 2 ; Hj^ P (U ; //)) is defined analogously. 

Definition 2.1. A function u E L p (t\, t 2 ; H,' p (U; /x)), where U is an open 
and bounded subset of M. d , is a weak solution of equation (jl.ip in U x (r%, r 2 ) 
if for any compact SI C U and for almost every n < t% < t 2 < t 2 it satisfies 
the integral equality 



[u q (x,t 2 )r](x,t 2 ) - u q {x,ti)r){x,tij\ dfi 

o (lv.p-v..v,-*£)*-o 



for every 77 G Cq°(S1 x (h,t 2 )). 
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Here the boundary terms are taken in the sense of limits 

rti+h r 

u(x, tx)r](x, ti) dfi = lim / / u(x,t)r)(x,t) dv 
n h->oj tl J n 

and ^ 

u(x,t2)r](x,t2) dft = lim / / u(x,t)rj(x,t) dv. 
< h ^ Jt 2 -hJu 
The following notion of weak solution given by Steklov averages is technically 
more convenient to work with since it involves the discrete time derivative of 
the solution and it is equivalent to the previous definition, [5, p. 101-102]. 

Definition 2.2. A weak solution of (jl.ip is a measurable function u E 
L foc(n, t% H]f c (U ; fi)) satisfying 

f l^P^T] + (\Vu\ p - 2 Vu) ■Vr j \dv = 0, (2.2) 
Jnx{t} { ot V J h ) 

for every compact O C U and for almost every t\ < t < T2 — h and for all 

Next, we recall a few definitions and results from analysis on metric mea- 
sure spaces. The measure \x is doubling if there is a universal constant 
Dq > 1 such that 

fx(B(x;2r))<D 0fJ ,(B(x,r)), 
for every B(x,2r) C £1. Here 

B(x, r) := jy G R d : \y - x\ < rj 

denotes the standard open ball in M. d with radius r and center x. 

The dimension related to the doubling measure is d^ = log 2 -Do- Observe 
that for the Lebesgue measure this is dc = d. Let < r < R < oo. A simple 
iteration of the doubling condition implies that 

SCC*^, (2.3) 



fi(B(x, r)) \ r y 

where C depends only on the doubling constant -Do- 

Recall that, for an a € (0,1], we say that the measure /i satisfies the 
a-annular decay property if there exists a constant c such that 

fj,(B(x, r) \ B(x, (1 - 5)r)) < c5 a fi{B(x, r)), (2.4) 

for all B{x,r) C £1 and for all 5 € (0,1). Obviously in our setting this 
property holds, this follows directly from the doubling property, in particular 
from (|2.3p . It was proven by Buckley [3] that any length space has the 
chain ball property and, for any metric space equipped with a doubling 
measure, the chain ball property implies the a-annular decay property for 
some a G (0, 1] (see also [1]). 

The measure is said to support a weak (q,p)-Poincare inequality if there 
exist constants Pq > and r > 1 such that 

/ „ \ Vp 



\u - u B ( x ,r)\ q dfx } <P r f |Vu| p dH , (2-5) 

B(x,r) I \JB(x,tt) J 
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for every u € H,'^ (0; /x) and B(x,rr) C £1. Here, we denote 

ub(x,v) = { ud[i= —j—r — rT / udfi. 
Jb(x,t) n(B(x,r)) J B (x,r) 

The word weak refers to the constant r, that may be strictly greater than 
one. In M. d with a doubling measure, the weak (q, p)-Poincare inequality with 
some r > 1 implies the (q,p)-Poincare inequality with r = 1, see Theorem 
3.4 in [12] . Hence, we may assume that r = 1. 

We will assume that the measure supports a weak (1, p)-Poincare inequal- 
ity. On the other hand, the weak (1, p)-Poincare inequality and the doubling 
condition imply a weak («;,p)-Sobolev-Poincare inequality with 



l<P<dn, 



d^-v ' ^ ^ ^ (2.6) 
,2f>, P>d^, 

where d^ is as above. For the proof, we refer to [12J. 

For Sobolev functions with zero boundary values, we have the following 
version of Sobolev's inequality. Suppose that u € HQ ,p (B(x,r); y). Then 

i/p 

\u\ K dy) <Cr\4- \Vu\ p dy\ . (2.7) 

'B(x,r) J \JB(x,r) J 



(£ 




For the proof we refer, for example, to [18]. Moreover, by a recent result 
in [16], the weak (l,p)-Poincare inequality and the doubling condition also 
imply the (1, g)-Poincare inequality for some q < p. Consequently, we also 
have the weak (ft, g)-Sobolev-Poincare inequality which implies, by Holder's 
inequality, the (q, g)-Poincare inequality for some q < p. For more informa- 
tion on Sobolev- type inequalities see [231 12"4] . 

Indeed, the Euclidian metric does not play a crucial role in our method 
and can be replaced with another metric provided that the space endowed 
with this metric has the annular decay property. As mentioned before, this 
property holds true for length spaces, i.e. metric spaces in which the distance 
between any pair of points is the infimum of the lengths of rectifiable paths 
joining them. Moreover, the fact that the theory of Sobolev spaces is based 
on the notion of weak derivatives brings about the need of an alternative 
way of defining weak gradients in a general metric measure space. However 
this can be overcome by employing the so called upper gradients. 

Our main result is the following theorem. It still holds for general equa- 
tions with a more general principal part A 

IguiW - div (A(x, t, v, Vv)) = 0, 
where the Caratheodory function A satisfies the structure conditions 

A(x,t,v,r l )-ri>A Q \r l \ p , (2.8) 
A(x,t,v,r ] )<A 1 \r ] \ p -\ (2.9) 

for almost every (x, t) £ W 1 x M and every (it, rf) € R x M n , for some constants 
Aq,Ai > 0. 



LOCAL HOLDER CONTINUITY FOR DNPE 



5 



Theorem 2.3. Let p > 2 and < q < 1. Assume that the measure [i 
is doubling and supports a weak (l,p)-Poincare inequality. Then any non- 
negative weak solution of equation (jl.ip is locally Holder continuous. 

We will work with parabolic cylinders built upon balls. Let (xo,to) be a 
point in the space-time domain. The cylinder of radius r > and height 
s > 0, with vertex at (xo,to), is defined as 

Qx ,t (s,r) := B(x ,r) x (t - s,t ). 

We write Q(s,r) to denote Qo,o{s,r). Moreover, we shall use the notation 
5Qx ,t (. s i r ) = Qx ,t {^s,Sr) and SB(xo,r) = B(x ,5r). 

In the sequel, we shall refer to data as the set of a priori constants p, q, 
d, Do, and P$. We will suppose that the weak solutions of equation (jl.ip 
are locally bounded and we will assume 

for some positive constant M. This assumption will be used without further 
comments hereafter. 

We will need the following algebraic lemma, the so called "lemma on fast 
geometric convergence" (cf. [5]) to establish our result. 

Lemma 2.4. Let (Y n ) n be a sequence of positive numbers satisfying 

Y n+1 < Cb n Y^ +a , 

where C, b > 1 and a > 0. Then (Y n ) n converges to zero as n — > oo provided 

Y < c^^b 1 -* 2 '. 

3. Constructing the geometric setting and fundamental 

estimates 

Let u be a nonnegative weak solution of equation (jl.ip in U x (n, T2) and 
let if be a compact subset of J7 x (t\,T2). The Holder continuity of u at a 
point (a^o^o) i n K follows via an iteration process applied in a sequence of 
nested and shrinking cylinders with vertex at that point. In each step of the 
iteration process, we prove that the oscillation of the solution reduces in a 
measurable way as we suitably decrease the size of the cylinder. Finally, the 
oscillation converges to zero as the cylinders shrink to the point. Because 
it is always possible to translate the equation, without loss of generality, 
we restrict the study to the origin. Moreover, since we focus on the local 
interior regularity it is enough for our purposes to assume K is the cylinder 

K := Q(2R 2 ,2R), R>0. 

For the initial cylinder K, let 

H~ < ess inf u and /x + > ess sup u, 
K K 

and define 

to := /i + — 
We choose pT small enough so that 
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holds. We may assume that uj > 0, because otherwise there is nothing to 
prove. We also assume that, without loss of generality, u < 1. 

For a suitably chosen sequence {Q 1 } of cylinders, we will construct a 
nondecreasing sequence and a nonincreasing sequence {fif} so that 

< essinfu and fif > ess sup u. (3-2) 

and 

Ui := fj,f - [i~ = a l u, i = 0,l,... 
for some a € (0, 1). As a consequence of (|3.2j) we have 

essoscu < cjj. 

The actual proof will proceed by induction. We will assume that in Q l the 
two-sided bound f|3.2[) holds. Then we shall built the cylinder Q t+l in such 
a way that (|3.2p is verified with i + 1. 
For 5 > 0, sufficiently small, let 

^ = 5*11, i = 0, 1, (3.3) 

We will work with a sequence {Q 1 } of nested and shrinking cylinders of the 
form 

g J : =Q( Ci <,^,), Ci = o;r 1 (^J , i = o,i,..., 

where A > 1 is a constant to be fixed later depending only on the data. 
In particular for the initial cylinder Q° = Q(cqR p ,R) we choose fj, as 

l4 = ^ (3.4) 

which immediately gives 



ujq = oj and cq = uo q 
We assume, without loss of generality, that 



i?<2-V + H. (3.5) 

Indeed, otherwise we would have co < 2~ X ^R P , with /3 = [1 + (1 - q)/(p - 
2)] _1 , in which case there is nothing to prove because then the oscillation is 
comparable to the radius. Assumption (|3.5p guarantees the inclusion 

Q(coR p ,R) C Q(R 2 ,2R) 

and thus implies the starting relation 

ess sup u <ojq := u. 
Q° 

In order to assure the inclusion Q l+1 C Q l for all i > 1 it is enough to 
assume 5 < atP-v-V/P. 
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3.1. Energy estimates. We introduce the auxiliary function 



j((u-k) ± )=± / (e /q -k) # 

Jki v 7 ± 

=± q r^-k)^- 1 ^ 

Jk 



An— k)± 

Jo 



for which 

d d(v q ) 

^J((u-k) ± ) = ±^-l( u -k) ± . (3.6) 
In what follows, we will need the following estimates. For p > 2, we have 

J((u-k)+) = q k+ (k + 9 ~^dC 



f (u— fc) + 

iqk"- 1 I (<% (3-7) 



and 



2 

(«-*0h 



J{{u-k) + )>qu q - 1 / £ci£ 
■> 

-l(u-*£ 



(3.8) 



> git 9 " 



2 

On the other hand, for the minus case we obtain the following estimates. 
J{(u—k)-) = q "'(k-tf-^dt 



,"(»-*)- . (3.9) 

< q(u-k)- / (Jfc 



and 



< k q (u-k)- 

J{(u-k)^)>qk q - 1 ^ £dt 
>qk q -^- 



(3.10) 



2 

We will deduce the fundamental energy estimate over the cylinders x 
(h,t 2 ) CUx (n,r 2 ). 

Lemma 3.1. Let u > be a weak solution of (jl.ip and let k > 0. Then 
there exists a constant C = C(p) > such that 

ess sup / J((u—k)±)(p p dfi + / / | V(u— k)±ip\ p dv 
t!<t<t 2 Jn Jti Jn 

< C / (u-kf ± \Vv\ p dv + C f 2 [ J((u-k) ± )^- 1 dv, 
Jti Jn Jti Jn \ at J + 

(3.11) 
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for every nonnegative if G Cq°(£1 x (ii, ^2))- 

Proof. We take 77 = ± — k)± (p p , where Uh = [{u q )h] 1 ^ q , as a test function 
in (j2.2|) . which is in the Sobolev space H ' p (p.) and only admissible after an 
approximation, and integrate in time over (ti,t), for t € (ii,^) to obtain 
i 



= 

ti n 



d((u q )h) . f lv? ,p-2 r 1 v- 

77+MVul Vit) ■ Vry 



The first term on the right hand side can be estimated from below by us- 
ing (|3.6p . integrating by parts and finally letting h — > as 



<H '"' )Q [±(u h - k) ± <f?] du= I (J ((u h - fc)±)) ^ d!/ 



dt 1^ 1 J J dt 

ha t\ q 

t 

'dip 
W 



j J((u-k)±)<pPdti-p j jj({u-k) ± )if- l (^ 



dv. 



nx{t> h n 

Concerning the second integral on the right hand side we first let h — >• and 
then use the estimate 

|V(« - £;) ± | p - 2 (±V(u - k)±) ■ V(±(u - k)±p p ) 

> |v(« - fc)±|V -p|v(« - fc) ± | p -V _1 (« - *0±|V¥>|. 

The last term is estimated further by Young's inequality as 
-p|V(u-fc)±p'-V~ 1 («-*0±|V¥'| 

> -||V(«-%|V - C(u-k) p ± \Vip\ p . 

Hence the estimation of the second integral on the right hand side reads as 

t 

\Vu\ p ~ 2 VuJ -Vrjdv 

h n 

t t 



\J J\V(u-k) ± <p\* dv-C( V )J J{u-kf ± \V V \ p dv. 

tain (fXTTD . 



> 

2 

ti si ti a 

Since t E (^1,^2) is arbitrary we can combine both estimates to ob- 



□ 



3.2. Logarithmic Estimates. We introduce the function 

faW := »(af, («-*)*,„) = (m ( c + H "_\ u _ k) J ) + - 

for which we obtain the following logarithmic integral inequalities. These 
estimates are technical tools used to forward information in time. The proof 
is somehow formal since it involves the time derivative of v. A rigorous 
justification requires the use of Steklov averages in the spirit of [5j p. 101- 
102] (see also [21]). 
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Lemma 3.2. Let u > be a weak solution of equation Then there 

exists a constant C = C(p,q) > such that, for p > 2 and < q < 1, we 
have 

fc 9-1 ess sup f ip 2 L(u)(x,t)(p p (x)dfx< [ u q ~ l ^ 2 _(u)(x, h)(fP(x) dfi 
ti<<<< 2 Jn Jn 

+ c f 2 f v^)l(voV)Mv<^ 

h\ Jn 



and 



esssup / ^ + (v)(x,t)ip p (x) d\x < / ip + {v)(x, t\)ifP{x) d\i 
t!<t<t 2 Jn Jn 



,-t 2 



(l-q)(p-l) 



+ C I I v^r-^ + (v)\& + ) (v)\ 2 -P\Vip\Pdv, 
Jt\ Jn 

where if € Cq°(Q,) is any time-independent test function and v = u q . In the 
plus case we assume that u is strictly away from zero. 

Proof. We start with the minus case. Choose 

in the definition of weak solution and integrate in time over (ti,t) for t E 
(tiftz). Observe that 

(^y = 2(i+v-)(v-) 2 . 



(3.12) 



The estimate of the parabolic term reads as 

* r du q , , r* f d 
—■q-[u) dv 



ti Jn 



h -Jn 9t j kq 



ri-ts 1 /*) ds dfj, 



n Jki 




1 1 



r/_(s)s q 1 dsdfjL 



n Jk 



(3.13) 



Integrating by parts we obtain 



n-(s)s q - 1 ds= / (^ 2 _(s))'s q ' 1 dsi P p 



</[^( s ) s <?-i]«_ ((? _i) j ^_( s )s q - 2 ds^ 

rk 



= iP 2 _(u)u q - 1 L P p -(l-q) / 7pl(s)s q - 2 dsip p . 

J u 

From the above equality we get 

r ] ^(s)s q - 1 ds > ^ 2 _{u)k q - l ^ p 
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and 



77_(s)a ?_1 da < ^(u)u« - V- 
Using these estimates in (|3. 13[) gives 



ft r g u i 



tp 2 _(u)(x,t)(p p (x) d/i 



n 

u q - 1 ^ 2 _(u)(x,t l )(p p (x)dfi, (3.14) 

n 



for all i G (ti,t2)> 

Concerning the remaining term, we first let /i — > and by using (|3,12|) 
together with Young's inequality, we obtain 

\Vu\ p - 2 Vu • Vt?_ = \Vu\ p - 2 Vu • v((^i(«))V) 

> 2|Vuj p (l + V - 2p|Vu| p ~V-V'-</~ 1 |V9?l 

> |Vu| p (^'_)V " C^_(^) 2 " P |V^| P , 

almost everywhere. The claim for the minus case follows from the previous 
estimate and (|3.14p . 

In the plus case we choose 

n + (v) = ^l( v ))^, 

where v = u q . Notice that (|3.12p continues to hold for this choice of test 
function. After integrating in time over (t\,t) for t € (ti,^)) the estimate 
of the parabolic term follows as 

ti Jn at 

h 

By using (|3.12|) together with Young's inequality, we obtain 
\Vu\ p - 2 Vu • Vr? + = \Vv l / q \ p ' 2 Vv 1 ' q ■ V((^(v))V p ) 



n 



qP- 1 



p (i-g)( P -i) „ 
+ ^ZT« ^ |Vf| p_2 Vw{2V'+V'V v ^} 

> -C7(p,g)t; ii ^ ii ^(i;)(^(«)) 2 - p |V^| p , 

almost everywhere, from which the claim follows. Note that here we also 
use the fact that f(x) = < q < 1 is Lipschitz for strictly positive 

values of x. Therefore by the chain rule we get 

i / 1 1-9 

Vv^ q = -v~ 

q 

since in this case we assume that u is strictly positive. □ 
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4. Continuity analysis 

Suppose that u > is a solution of (jl.ip in Q % = Q (ciR^,Ri) and we 
have 

lij < essinf u < ess sup u < [if, fif — =: Wj. (4-1) 
Q (> 

Assume that [if satisfy 

4m," <"i = t4~lh' ( 4 - 2 ) 
The case where this assumption does not hold will be treated in the Sub- 
section 14.31 at the end of this section. 

Inside Q 1 , we consider sub cylinders of smaller size 

which are contained in Q % whenever 

_ 2 A(p-2)j g_ < f < q ^ 

wf 9 

Since A can always be arranged so that iVi = Cj/dj = 2( A ~ 1 )( p_2 ) to be an 
integer Q l can be viewed as the union, up to a set of measure zero, of N-i 
disjoint cylinders each congruent to Qo,t* (diR p ,R). 

The reduction of the oscillation is based on the analysis of an alternative 
(see [3 [26]). For a constant ao € (0, 1), that will be determined depending 
only on the data, either 

The First Alternative. 

there is a cylinder of the type Qt*(diR?, R{) 

u({(x,t) € Q t *(diR%,Ri) : u(x,t) < #r + Wi /2}) 



or, since u T — -f = /j, + -f 



v(Q t * (d i R p l ,R l )) 

2 



< a , (4.4) 



The Second Alternative. 

for every cylinder Q t *(diR p , Ri) C Q l , we have 

f({(aj,t) € Q t *(diR?,Ri) : u(x,t) < /i+ -^/2}) 



< 1 - a . (4.5) 



i/(Q r (diR?,Ri)) 

In both cases, we will prove that the essential oscillation of u within a smaller 
cylinder decreases in a measurable way, measurable so that we can derive 
a modulus of continuity. The constant ao will be fixed in the course of the 
proof of Lemma 14.11 

4.1. The First Alternative. Now we assume that the first alternative 
holds. 

Lemma 4.1. There exists a constant ao £ (0, 1), depending only on the 
data, such that if (|4.4p holds for some t* as in (j4.3j) then 

u(x, t) > + ^ in ^Q t * (diR 1 ?, Ri) . 
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Proof. Define 

R n = Y + ^TT' Qn = B n xT n = B(R n ) x (t* - diRP, t* 

and 

u i ( 1 

kn = * + T\ ¥ 

for n = 0, 1, ... . Choose cutoff functions (p n € C°°(Q n ), vanishing on the 
parabolic boundary of Q n , and such that < ip n < 1, (p n = 1 in Q n +i, 

- 

By assumption (|4.2p 

j < < Uj. (4.7) 

Here the first inequality is trivial and follows from the definition of k n . 

By Holder's inequality and Sobolev's inequality (|2.7p together with the 
doubling property of the measure, we have 



J Qn + l 



(u-k n )f- p/K)+p du 



< v(Q n ) 



.1 ( u -k n ) 2 t 1 ~ p/K)+p ^ 1 -P/^ du 



>(Qn+l) 

<c4 4 (u-fc n )|/ n ^ 4 ((u-k n ) ± <p n ) K dfi) dt 

J T n \J B n J \J B n J 

<CR P ( esssup-f (it-A^)±^d^tj + |V((it-fc n )±yn)[ p dz/, 

(4.8) 



where k is the Sobolev exponent as in (|2.6|) . Note that this estimate is 
independent of the choice of k n and Q n . 

Taking into account (|3.9p and (|3.10p . energy estimate (|3.ip in the minus 
case takes yields 

k^ 1 esssup + (u— k„) 2 _^ d\i + d{R p A \ V(u— k n )_(p n \ p dv 

Tfi J J 



13 n 



p-l ( dl Pn 



< C diRpj \u-k n ) p _ \Vvn\ p du + CdiRlkljiu-K).^-' dv. 

Qn Qn 

Using (|4.6|) and (|4.7p together with the fact that (u— k n )_ < almost 
everywhere the above inequality can be rewritten as 

ujf~ esssup -j-(u—k n )^_ ifP n d[i + diR^-f |V(ti— k n )-ip n \ p dv 
T n J ' J 

B n Qn 

<c2^dAm p +^ l ^ v{An) 



2 di J v{Q n ) ' 
where 

A n = {(x, t) EQ n : u(x, t) < k n } . 
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Thus, by using the definition of di, we conclude with 

ess sup/ {u-KUld^ < CV* m 2 ^ (4.9) 



and 



K j M u- kn) _ vTdv <C2^(^Y^ y (4.10) 



Qn 

The last step is to estimate the left hand side of (|4.8p from below 



JO. 



(u-k n ) 2 ^- p/K]+p du 

n + l 

>f (« - k n ?l l ~ plK)+P X{{u-k n+1 )->v} dv 

J Qn + l 

> \k n+1 - K^-v'^ I x {(u - fcn+l) ->o} du 

jQn + l 



(4.H) 



> 



Ui \2(l-p/«)+p v(A n+1 ) 



,2™+ 3 / u(Q n+i y 
Inserting (|4U|) . (|4"?T0"|) and (|4~TT|) into (pL"8j) we conclude 

VAAj}±A < ^ f 4^4 T P ' K , n = 0,1,2,.... (4.12) 



By setting 

i/(A 



Y, 



u{Q n y 

we obtain the recursive relation 

v <r r <n + 1 v 2 ~p/ K 

I n+l ^ > 

for some constant C, depending only on the data. We conclude, using 
Lemma 12.41 that if 



Y < c-Va-pAO+Mi-pAO 2 .-_ 



a 



holds then Y n — > 0, as n — > oo. But since this condition is precisely the 
assumption of the first alternative ()4.4p for the indicated choice of cto the 
result follows. □ 

Now we impose the logarithmic estimates and make use of the previous 
lemma to prove that the set where u is close to its infimum can be made 
arbitrarily small within the time interval —9 < t < 0, where 

-O^f-di^J (4.13) 

for some t* as in (14.31). 



Lemma 4.2. Assume (|4.4p holds for some t* satisfying (|4.3p and let 9 be 

given as in (|4.13p . For every number ot\ S (0,1), there exists s* G N, 
depending only on the data, such that 

n(lxE^B{Ri) : u ( x ,t)<K+^X\ <a xl ,{^B{BM , Vie (-9,0). 



14 
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Proof. Our aim is to apply Lemma 13.21 to the function (u — with 



k = fi i + 



UJ; 



and H,, 



4 ' 2™+ 2 ~~ K ! 

in the cylinder (9B%, R4). Here n G N will be determined later depending 
on the data. Observe that, by Lemma |4.H u(x, —9) > k in ^B (R{). Thus 



^(u)(x,-9) = 0, xE-B(Ri). 



On the other hand, since (u — k)-. < we get 



UJ-: 



and 



2-P 



n In 2 



(^T-(--^)- + C r 2 <(y 



Choose </3 G Cq° (2-B (-Rj)) independent of time and with properties < 92 < 
1, = 1 € |J3 (J?i) and 

|V¥>| < 



In the set 



R, 



xe-B (Ri) :u< N + — 



ip-{u) can be bounded from below as 
ip-(u) > In 



2 -(n+2) + 2 -(n+2) 



> (n - l)ln2. 



Note that (02]) implies fe < ^ . 

Now we apply Lemma 13.21 in the minus case to the function (it — and 
use the above estimates to obtain 



^- 1 (n-l) 2 (ln2)V(<jxG^(i? i ) : « < ^ + ^ 



1 



< k q - 1 



< c 



ess sup 

-8<t<0 J ±B(Ri) 









(V>-)'(«) 







2-P 



<Cn(ln2)-^0/z (^R(Ri 

<Cn(ln2)(g) 2 ^ W rV (^(^))- 

In the last inequality we used the fact < CjiZf . A simplification in the 
above inequality gives 



,i[{xe^B (Ri) :u<Hi + y 



n 



To conclude, we choose s* = n + 2 with n > 1 + ^-2 A( - P 2 ). 



□ 
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The information of this lemma will be used to show that the conclusion 
of Lemma 14.11 holds in a full cylinder that includes the origin. 

Lemma 4.3. Assume (|4.4p holds for some t* satisfying (|4.3p and let 9 be 

given as in (|4.13p . There exists si £ N, depending only on the data, such 
that 



Proof. Let 



and 



Rn — "^" + ^^3 > Qn — Q(9,R n ) 



, _ CJj / 1 

L = M,- H tt 1 H 

" ^ 2 Sl + 1 V 2" 



for n = 0,1,.... Take cutoff functions <p n (x) € C°°(B n ), where B n = 
B(R n ), vanishing on the parabolic boundary of B n , and such that < ip n < 
1, (p n = 1 in S n+1 , 

C2 n 

|V<y? n | < 



Observe that k n < u>i by assumption (|4.2|) and u(x, —0) > fJ>i + ^ > k n in 
-B„ C ^B(Ri) by Lemma |4 . 1 1 which implies —0) = in £? n . Taking into 
account these estimates, (|3.9p and (|3.10p . energy estimate (|3.ip in the minus 
case takes form in this setting as 



5-1 

uS\ ess 

-0<t<O 



Qn 

Using the facts 



sup -f (u— kn) 2 _ifP n d\x + 9 4 \V(u—k n )-f\ p dv 
a<oJ J 

<C9j(u-k n ) p _\Vy n \ p dv. (4.14) 



(u — k n )- < — - and 9 < c;R p 
~ 2 S1 ~ 



together with the estimate on |V<^ n | we arrive at 

ess sup l(u-k n f_^ < CI** (^)V- Sl)(p " 2) ^4 
-e<t<oJ v 2 s i/ v{Qn) 

and 

Qn 

where A n = {(x,t) € Q n : u(x,t) < k n }. Observe that the estimate (|4.8p 
continues to hold in the setting of this lemma and the left hand side is 
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estimated from below as 

(u-k n f^- p/K)+p du 



£ 

Jo 



n + l 

>f ( u ~ k n ?} l ~ P,K)+P X{(u-k n+1 )->o} dv 

Qn+1 (4.15) 
> \k n+1 - k n \ 2 ^^ +p / X{(u-k n+1 y> 0} dv 



> 



J Vn + l 

u . v2(l-p/«)+P u(A n+1 ) 



W*"+*J u(Q n+i y 

Substituting these last three estimates in (|4.8p and assuming that s± > A, 
we obtain 



\2 



u(Q n +x) \v(Qn) 

We set 

as before and rephrase the above inequality as 

v <r r <n + 1 v 2 ~p/ K 

J-n+l _i ^ I n 

Define 

a x = c-i/(i-p/«)+i-(i-p/«) 2 . 

We conclude, using Lemma 12.41 that if Yq < a± holds then Y n — > 0, as 
n — > oo. Now we apply Lemma 14.21 with such an a\ and conclude there 
exists s* =: s\, depending only on the data, such that 

n(LejB(IU) ■■ u{x,t) < fir + < axfi Qfl (rA , Vt G (-0,0), 

which is exactly Yo < ao- The result is concluded since Y n — > 0, as n — > oo 
implies that A n — >• 0, as n — >• oo. □ 

We finally reach the reduction of the oscillation of the solution for the 
first alternative. 

Corollary 4.4. Suppose that ([321) and hold in Q i = Q(aR P ,Ri). 

Assume also that (|4.4p is verified. Then there is a constant aj € (0,1), 
depending only on the data, such that 

essosc u < cjWj. (4-16) 



Proof. Observe that 



,/,l§y <e = -t* + d i (^ P 



since t* < 0. Therefore 

ess inf u > ess inf u > a- H — - . 

iO(*fl?A) 2S1+1 

Put ^ 

Mw-i = and A»i+i = MT + ^ITr 
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Then we have 



1 



and the corollary follows with 07 = 1 — 2S1 1 +1 G (3/4, 1). □ 



4.2. The Second Alternative. Now we analyze the second alternative. 
Assume (|4.5p holds for all cylinders of the type Qt* (diR?,Ri), where t* is 
as in (B~3l). 



Fix a cylinder Q t * (diR?,Ri). We deduce from (|4.5p that there exists a 
time level t° G (i* - d^f, t* - ^dj-Rf ) such that 

M({s €5(120: n(x,t°)>^+-^/2}) < f /j"^ ) KW)- (4-17) 

In fact, if g37D is violated for all (t* - djUf ,f - ^di-flf), we would get 
f({(M) G Q t * (diR p ,Ri) : «(x,t) > fj+ - Ui/2}) 



> 



t*-(a /2)d t RP 

Li{{x£B(Ri) : «(s,t) >$-Ui/2}) dt 



> (1-aoMQ (d i R 1 i,R l )), 

which contradicts (|4.5p . 

The next lemma asserts that a similar property still holds for all time 
levels in an interval up to the origin. 



Lemma 4.5. Assume (|4.5|) is verified. There exists s, depending only upon 
the data, such that 

li{{xeB{Ri) : u(x,t) >i4-2~ s u H }) < (l - ^ ^B(Ri)), 

for almost every t G (— (ci/2)R?, 0). 

Proof. Let 

c = 



>s+l 



and set Q := B{Ri) x (t°,i*). We shall use Lemma 1331 to forward the 
information in time. Set v = u q and recall the definition 



il>+(y) = V{H+,{v-k) + ,c) = ln + (- 



c + H+ -(v-k). 
From (v — k) + < (iOi/2) q = Ht we obtain 



^ + (v)<\n-^- = sq\n2. 



We continue with the estimate 

W+(v)?~ P < i H k+ c Y~ 2 < 2^- 2 ^f p ~ 2) . 
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Denote B := B{Ri) and let 93 £ Cq°(B) be a time-independent cutoff func 

c 



tion and has the properties < ip < 1, (p = 1 in (1 — <5)£> and |Vc/?| < jS-, 



where < <5 < 1 is to be determined later. Define 

s={xe(i-s)B: v>(,tr-(^y} 

and 

tf = {*€(l-d)B : v>(4-^) q } 

and observe that S' C S" for g € (0, 1). 
In the set 5 we get 

ib+(v) > In 1 2 ; . ^ > In ) 2 ' a > (sq- l)ln2. 

Now we apply Lemma 13.21 in the plus case with these choices to conclude 

(sq- l) 2 (ln2)V(5') < ess sup / ip$(v)(x,t)(fp(x)dn 

t°<t<t* Jb 

< I ^ + (v)(x,t°)^(x)d t i 
Jb 

+ C f [ v^^^ + {v)\^ + )'{v)\ 2 -v\V^dv 
Jt° Jb 

<( Sg ) 2 (ln2) 2 -i^M£) 

1 - Qo/2 

for almost every t £ (t°,t*), where C depends only upon the data. Observe 
that in the third inequality we have used (|4.17p and 

V = U q < (2L0i) q , 

which follows from (14.21). We use \t* — t°\ < djFtf to obtain 



By the annular decay property (|2.4|) . we have 
fi({x € B : v(x,t) > (fif - 2-( s+1 )^) 9 }) 

< fi(B \ (1 - 8)B) + fx({x e (1 - S)B) : v(x,t)> (fif -2-( s+1 ^i) q }) 

< C5 a f i(B) + n(S'), 

for almost every te(t°,t*). For the first term, we choose 5 small enough so 
that 

24 

and for the second term we use the previous estimate. By choosing s large 
enough so that 

1 - «0 {sq) 2 < _ Qo 
1 - a Q /2 (sq -I) 2 ~ 3 
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and 



C{sq) a 



6P(sq - l) 2 ~ 24 
hold, we get 

H ({* € B : v(x,t) > - 2^ +1 )u; i )''}) < 1 - ^ (4.18) 
for almost every 

Since ()4.17|) holds for all cylinders of type Qt* (diR p , Ri) the conclusion (|4.18p 
holds for all time levels 

t>-(ci-di)B$- ^-diRl 

By choosing 

2 (A-l)(p-2) > 2 ^ 419 -) 

we get Cj/di > 2 — ao, which implies 



□ 



The information of the previous lemma will be used to prove a critical 
estimate for the reduction of the oscillation in the second alternative. 

Lemma 4.6. For every a.i 6 (0, 1) there exists m > s, depending only on 
the data and a.<i, such that 

v ({(x,t) £ Q (^RP,Rj) : u(x,t)> y+ -2- m Ui}) < 

v{Q{%rP,dir)) -° 2 ' 

Proof. Let 

E n (t) = {xe B(Ri) : u{x,t) > n+ - 2~ n Ui} 

and 

E n = {(x,t)eQ(^R P ,R^j : u(x, t) > ^ - 2~ n Ui}. 

Denote 

h = nt- 

and 

k = fif — 2~ n uj u 

where n> s will be chosen large. Here s is as in Lemma 14.51 Let 

h — k, u > h, 
iv = < u — k, k < u < h, 
0, u < k. 

We shall use the notation B = B{Ri). By Lemma 14.51 and the fact that 
n> s, for almost every t £ (— (ci/2)R p , 0), we have 

Li(xeB : w(x,t) = 0}) = € £ : u(sc,i) < fc}) > ^j^{B). 

Thus, we obtain 

w B (t) = -[ wdfi< (l-^){h-k) 

JBx{t} V 4 / 
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and, consequently, 

h-k-w B (t) > !j(h-k), 

for almost every t £ (— (ci/2)R p , 0). Using the (g, g)-Poincare inequality for 
some q < p (see (|2.5p and the subsequent remarks), yields 



(^-fc)V(^ + iW)<f-V / 

<CR q [ \Vw\ q dn = CR q t [ \Vu\ q dfi, 

JbxH} J E n (t)\E n +i (t) 



\w — WB(t)\ 9 dfj, 
'Bx{t} 



'Bx{t} J E n (t)\E n+1 (t) 

for almost every t € (— (ci/2)R P , 0). The constant (4/«o) 9 above is ab- 
sorbed into the constant C since it is already fixed in terms of the data in 
Lemma 14.11 Now we integrate the above inequality over time to get 

(h - k) q v(E n+1 ) < CR\ j \Vv\ q dv. 

J E n \E n +i 

Introduce a cutoff function tp 6 (Q (ciR p , 2R.{)) such that < tp < 1, tp = 
1 in Q ((ci/2)RP, Ri), if vanishes on the parabolic boundary of Q (ciR P , 2R{), 
and 

|V„|<£ and (%) 



Ri ' \dt) + - Ci R\ 
Now, Holder's inequality gives 

(h - kyu{E n+1 ) 

<CR q ( [ \Vu\p dv) v{E n \ E n+l ) l -i/v 

\JE n \E n+1 J 

<c(r? [ |V(u-fc)+|Vd^ v{E n \E n+1 ) l - q lv. 

Assume A > n(> s). The first factor on the right hand side can be estimated 
by Lemma \3. II in the positive 



|V(u-fc)+IV dv 



»-*0+ ( If ) (h ' 



<C f {u-k) p + \V^\ p dv + Ck q - 1 f 

jQ(c,R p t ,2R z ) Jq( Ci RP,2R 

In the second inequality we have used (u — k) + < p- < || and k > In 
the third one the doubling property of the measure v is used together with 
the assumption A > n. We obtain 



(^yVn + l) < C [„ (Q 



\1-<j/p 
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Finally, summing n over s, . . . , m — 1 gives 
(m - s)v{E m f/^ <c[v(q (|U?, J2< 



p/(.p-q) 



and hence 



c 



(m - s )(p-9)/p 
Choosing m large enough finishes the proof. 



□ 



Now we prove the main lemma of this alternative. We also determine the 
size of the cylinder Q l by fixing A and consequently q. 

Lemma 4.7. Assume (j4.5|) holds. Then the choice of X can be made so that 



u(x, t) < ill 



a.e. in -Q ( — , Ri 
2 V 2 



Proof. Define 



2 2 n+1 ' 



2 A+1' 



B n = B{R V 



and 



k n — 



2 A+1 



1 + 



1 



for n = 0, 1, . . . . Consider the cutoff functions (p n {%) S C^°(Q n ) vanishing 
on the parabolic boundary of Q n and such that < tp n < 1, tp n = 1 in <5n+i, 



C2" 

|Vc/?„| < — - and 
tu. 



dt 



< 



C2 pn 



(4.20) 



Observe that u < fxf = \i~ + < 2wj by assumption (j4.2j) . We also have 
the estimate for the levels k n > 2~ 1 uii. Allowing these estimates together 
with (j3.7p and (j3.8f) . the energy estimate ()3.ip in the plus case takes form 
in this setup as 



UJ 



1 esssup -j-iu—kr^if^d^L H — TrRnT |V(u— & n )_(^| p 



1 B, 



< c 



(u-fc„)_<^ 



2 ,.J>-1 



5f 



(4.21) 
(4.22) 



Using ()4.20p and the estimate (u — k n ) + < 2 x uj{ we arrive at 

Wj\2 v (A n ) 



sup I{u-k n f_^ n dn<C2^(^ 



and 



OSS 



A B 



RP n j \V(u-k n ). 



ip\ p dv < C2 n P 



2 X J u{Q n ) 
0Ji\P u(A n ) 



2 V u(Q n ) 
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where A n = {(x,t) € Q n : u(x,t) > k n }. The left hand side of (|4.8p . which 
continues to holds in the setting of this lemma as well, is estimated from 
below as 



I 

Jo 



(u-k n )^- p/K)+p dv 

1+1 

>f ( u ~ M- (1 ~ P/ '" )+P X{(n-A ;n+1 )_>0} dv 

Qn+1 (4.23) 
> \k n+1 - M 2(1 ~ p/K)+P /_ X {(u - fcn+l) ->o} dv 

> 



J Vn+l 

u . x 2(l-p/«)+p i/(^ n+1 ) 



^2 A +"+V ^(Qn+l)' 

Substituting these last three estimates in (j4.8|) and assuming that si > A, 
we obtain 

2-p/K 



KQn+l) ~~ \v(Qr 



We set 

y = v{A n ) 
v(QnY 

as usual and we find 

v <r r" l + 1 y 2_ p/' : 

Defining 

a2 = cr-i/ci-p/Kj+i-a-p/*) 2 

we conclude, using Lemma 12,41 that if Yq < a\ holds then Y n — > 0, as 
n — > oo. Now we set 

A = max{l H — , m}. (4.24) 

p — 2 

The choice trivially satisfies (|4.19f) and by the previous lemma we conclude 
that A can be chosen large enough so that is as small as we please. 

Consequently we reach the final result from A n — > as n — > 0. □ 

We obtain the reduction of the oscillation of the solution as a corollary of 
the previous lemma for the second alternative. 

Corollary 4.8. Assume that holds in Q i = Q(aR p ,Ri) and that fljgj) 
is satisfied. Assume also that (|4.5p holds. Then there exists a constant 
an £ (3/4,1), depending only on the data, such that 

essosc u < crnuji. (4.25) 

±Q(^?A) 

Proof. By the result of the previous lemma we have 

ess sup u < p, 

Set 
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Then we have 

essosc u = uf,-, — uZi = I 1 rrr I w< 

and the corollary follows with ajj = 1 — ^att € (3/4, 1). □ 

4.3. The case where (|4.2p fails. If assumption (|4.2|) does not hold, which 
can happen only for an index i > 1 due to the assumptions (|3.1|) and (|3.4|) . 
then we have 

4 M r >oo iQ <=> /i+ <5/x~ (4.26) 
This elliptic Harnack estimate implies that if this condition holds for an 
index Iq then it continues to hold for all indices i > iq. By means of the 
(|4.26p we obtain (4^^) 9_1 < ^f™ 1 and this implies 



We redefine i?j as 
Therefore, we have 



- 4 (i_,)/p ^ U iv 



In this case we work with the following scalings factors 

-a = Kr 1 (It) 2 " and < = W" 1 (t ) 2 " p • 

The only significant change appears in the proof of Lemma 14.11 Instead of 
(|3.9p we use 

j((u-*0_) = q [ (u 

Jo 

r(u—k)— 

<qu q ~ l / idi 



< qu q 



2 

Observe that this estimate can be used in this case because thanks to the 
assumption (|4.26p we have u > > i > io, which implies that solution 
u is strictly away from zero. 

Now the whole process can essentially be repeated in the same way, the 
constants popping up along the calculations may change but they continue 
to depend only on the data. 

We finally prove the Holder continuity of u. Theorem l2.3l is an immediate 
consequence of the following theorem. 

Theorem 4.9. Suppose that u is a nonnegative weak solution of equa- 
tion (jl.ip in Q x t(R 2 ,R). Then there are positive constants C and a, both 
depending only on the data, such that 

( Q \ a 

ess osc u < C I — I ess sup u 

Qx,t(e p ,e) Q x , t (R 2 ,R) 



2d 
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for all < g < R. 

Proof. After translation, we may assume that (x,t) = (0,0). Take a := 
max{fj/, a n } and 5 := (T (p-<?-i)/p2-( 3 +( a - 1 )(p- 2 )/p). Observe that this choice 
of 8 gives 

C +1 (^< rf ,(f) P <|(f)". 
Then, by Corollaries 14.41 and 14.81 we have 

essosc u < Wi+i := <r 4+1 a;o = <r 4+1 w, i = 0, 1, . . . . 

Q(c I+ i(5»+ii?)P,5»+iR) 

But since Cj > 1 for all i = 0, 1, . . . , recall the assumption uj < 1 we finally 
obtain 

ess osc u < Ccr l uj, i = 0, 1, 2, . . . . 

Q((5 i i?)f,5 i _R) 

From this, the result follows in a standard way (cf. 0(26]). □ 
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